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Abstract 

High performance aircraft of the future will be designed lighter, more maneuverable, and operate over an ever 
expanding flight envelope. One of the largest differences from the flight control perspective between current and 
future advanced aircraft is elasticity. Over the last decade, dynamic inversion methodology has gained considerable 
popularity in application to highly maneuverable fighter aircraft, which were treated as rigid vehicles. 

This paper is an initial attempt to establish global stability results for dynamic inversion methodology as applied 
to a large, flexible aircraft. This work builds on a previous result for rigid fighter aircraft and adds a new level of 
complexity that is the flexible aircraft dynamics, which cannot be ignored even in the most basic flight control. The 
results arise from observations of the control laws designed for a new generation of the High-Speed Civil Transport 
aircraft. 

Introduction 

The advanced aircraft under consideration is a next generation supersonic transport aircraft, which due to 
aerodynamic factors will be long and slender and because of economics must be as light as possible. These factors 
contribute to making this aircraft very flexible with the first few elastic modes lying within the bandwidth of the 
traditional flight control system. 

Over the last decade, dynamic inversion methodology has gained considerable popularity in application to 
highly maneuverable fighter aircraft' which makes it a candidate that might benefit the advanced highly flexible 
aircraft of the future. This methodology has been modified and successfully applied to a new generation High- 
Speed Civil Transport vehicle 2,3 . However, there is very little applicable theory available on the stability of dynamic 
inversion based closed loop systems. Some recent work has been done by Morton, et. al. 4 focusing on the global 
stability of a rigid fighter aircraft. 

This paper is an initial attempt to establish stability results for dynamic inversion methodology as applied to a 
large, flexible vehicle. The method of reasoning follows the work of Morton, et. al.. This work builds on the results 
in reference 4 and adds a new level of complexity that is the flexible aircraft dynamics, which cannot be ignored 
even in the most basic flight control. The results arise from observations of the control laws designed for a new 
generation of the High-Speed Civil Transport aircraft that were described in some detail in references 2 and 3. 

In order to get initial results, the problem has been simplified for analytical work while retaining the essential 
characteristics. The simplification involved considering longitudinal dynamics with a single elastic mode and a 
control law designed for flight control only. In addition, throughout this paper the analysis considers the inner loop 
of the dynamic inversion only, i.e., the y*' toy; portion. It is important to note that the nature of y d ‘ 5 impacts the 

overall closed loop dynamics but will not be discussed here. 

The paper is organized as follows. The first section describes the equations of motion for this class of vehicles, 
followed by a definition of equilibrium set in section two. Section three presents the dynamic inversion results. 
Section four discusses stability results obtained to date, followed by concluding remarks. 


System Equations of Motion 

The equations of motion for an HSCT class vehicles (illustrated in figure 1) are rather involved (see reference 
3). They contain the usual aerodynamic forces and moments acting on the so called rigid vehicle and the actuator 
dynamics. In addition, equations for flexible body dynamics, their interaction with the traditional aerodynamics, the 
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interactions among themselves, and influence on the control surface activity are all also included. In this paper we 
limit our discussion to standard rigid body dynamics, one flexible mode, and interaction between the two. In 
addition, the control surfaces have been combined into one variable. Hence, a somewhat simplified longitudinal 
equations of motion for an HSCT type aircraft are given below 
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where 

m = vehicle mass, m n = generalized mass associated with an elastic mode 

GW * = mass distibution, T = thrust, M = Mach number, V = speed 

c - mean aerodynamic cord, S = planform area, q ~ dynamic pressure, p = density 

This paper considers a system with two control effectors, which are a traditionally low-bandwidth thrust and a 
relatively much faster acting elevator control. Though not explicitly written as a control, mass distribution GW *, 
can certainly be considered one. The GW * can be changed by shifting fuel around the vehicle as is done in the 
Concorde. In advanced version of a HSCT, there are no plans to employ it as such but GW* remains a degree of 
freedom in the system. 

The initent is to explore the stability of aircraft undergoing rapid maneuvering, which implies that the behavior 
of the fast states and parameters affecting them is of primary concern. This allows the treatment of thrust as well as 
the parametric dependence on GW* to be considered constant. In other words, the problem is reduced to a SISO 
control problem. Furthermore, in the subsonic regime we can legitimately drop the coefficient dependence on the 
Mach number. 

Equilibrium Set 

The system equations have the form 

x = f(x,#) 

where jc is a vector in R n and 8 is a vector in R m . Let U denote a set of allowed control values in R™ . Define the 
equilibrium set 

= {(*,*) | /(*,*) = o, seu}. 

Projecting M onto the first factor x we obtain A/, the set of equilibrium states 
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M={x\f(x,8) = 0, SeU) 

Note that M and M depend on the specified control limits. 

The equilibrium set of a system with m inputs is typically m-dimensional. The system modeled by (1) has three 
inputs ( T,GW*,S) and its equilibrium set M is three-dimensional. We choose to use T, 0 y and GW* as the 
coordinates on M Note that these are a mixture of controls and inputs. This choice allows us to eliminate the 
elevator control S/ from the equilibrium equations and leave them in terms of states that will later be shown to 
comprise internal dynamics of the closed loop dynamic inversion. In addition, choosing GW* allows for a unique r| 
for a given flight condition. 


Suppose that the triplet (T,6,GW*) is fixed, then at equilibrium rates and accelerations are zero giving us 
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Solving the pitching moment equation for the control value gives the following relation 
C M (a) + C Mn (a,GWy 7 
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Substituting for the controls in (2) and simplifying gives 
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The pitch rate dynamics disappear by design, which leave velocity and elastic mode dynamics. In the equilibrium, 
the force due to velocity acting on the elastic dynamics is balanced by the structural and aerodynamic forces and 
moments. Consequently, the dynamics of (u, w, if) define the rest of the equilibrium. 


The dynamics on the u-w-rj surface are 
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The forces of thrust and gravity must be balanced by the aerodynamic forces and moments. This formulation allows 
for a natural separation of dynamics due to rigid body and elastic interactions. 

We introduce the idea of equilibrium aerodynamic functions that will have direct connection to the aerodynamic 
forces required to maintain equilibrium state. 

Definition: For a system represented by equation (5), the equilibrium aerodynamic functions due to rigid body 
effects are C x (a), C z (a),and C n (a,V) defined by 

f ^ \ 


f C\(a) ' 

CM) 

v C 7 (a,^> 


G -c. 


C z -C, rfI 


c 


M 


c, 


c M 


c 


(a) 

(a) 


\f, ( i i 


V ( cos a + C 7 „ sin a ) - C, 


r?<V| 


Q,(«) 


the equilibrium aerodynamic functions due to the elastic effects are Cf(a ) , (\ ri (a) and C™(a) defined by 
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Hence the equilibrium aerodynamic force vector in the u-w-rj space at any equilibrium state is 
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Under practical circumstances an aircraft has a given flight envelope outside of which stability cannot be 
guaranteed or expected. This suggests dividing the analysis into two cases. 

1 . Assume unlimited control authority and work globally, or 
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2. Restrict analysis to a subset of states with adequate control authority. 


Due to space limitations, in this paper we proceed with analysis for the global case only. The restriction on 
control authority restricts the results to a subset that is discussed elsewhere. 

The system represented by (1) has three degrees of freedom and its three-dimensional equilibrium set M is 
characterized by solutions to the following equations 
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Remark: For a given triplet ( T,0,GW *), the equilibrium of system (1) is formed by all points 

(u,w,q,0irj,fj) = (w, h\0,<9, t/,0) 

where 

u- Vcos( a), w= Vs inf a), 


and (V, a,rj) is a solution of (6) corresponding to the given (T, 0,GW*). Thus given a triplet of (T, O y GW *), if (6) has a 
unique solution, then the system represented by (1) has a unique equilibrium point 


Dynamic Inversion 

The control problem can be stated as follows 


Problem statement: Given an equilibrium state x , determine a controller u - K(x,x) so that x is a global 
attractor for the system 

x = /(*, A'(x\x)) (7) 

Any global attractor must be an equilibrium state. Using dynamic inversion this problem is addressed for vehicle 
models having a unique equilibrium point for appropriately chosen engine thrust T and mass distribution GW*. 

The approach is to invert the rotational dynamics to a stable set of desired dynamics. Since the throttle is 
typically a low-bandwidth control that is not changed during dynamic maneuvers, it is left fixed throughout the 
construction of the controller K and the analysis. 

The philosophy behind the control law is for the vehicle to follow the pitch rate commands from the pilot and 
for the purpose of our analysis ignore the direct control of elastic modes. The structure of the dynamic inversion 
controller K is given by the following expression 

=K\(q cml -q) + K2{6^-0) (8) 


The desired dynamics are realized if the control surfaces conform to the following expression in the closed loop 
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Substituting the expression for controls (9) into the system equations (1) results in a closed loop system that readily 
separates into the following two subsystems. The controlled variable q and its associated equation gives 
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These equations are linear in nature and decoupled from the u-w- rj dynamics. The remaining dynamics of the u-w- rj 
subsystem result in 
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where the equilibrium aerodynamic functions have been substituted to simplify the expression. 

So proving the stability of the u-w- rj subsystem would prove the stability of the entire closed loop system (10)-(1 1), 
since the commanded variable is stable by design. 


Stability Result 

We begin by showing the stability of the internal dynamics system (11), followed by establishing that the only 
possible orbits of the residualized system, derived later in the section , are equilibria. Then we establish that the 
residualized system has a global attractor, and finally put all of these together to show stability of the closed loop 
system dynamics (10)-(1 1). 


Lemma 1: Assume the total aerodynamic drag, C n (a) = C D (a)-C DS (a ) — ~ (« ) > is always positive. Then there 
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exists a finite neighborhood D centered at the origin of the u-w-r/ surface into which all trajectories enter and 
remain. And the dynamic system described by (1 1) is asymptotically stable. 


Proof: Consider Lyapunov stability and the function V 2 + if as a candidate Lyapunov function. If the candidate 
function proves to be a Lyapunov function then we can make conclusions about the asymptotic stability of the 
internal dynamics given by (1 1). Since the function is a quadratic, it is easy to verify that it will satisfy the first 
criteria of a Lyapunov function. 

To satisfy the second criteria, from (1 1) we can compute the time rate of change of V 2 + i]" = u 2 + w 2 + z) 2 : 
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the V 3 term will dominate the candidate Lyapunov function. Hence, if the coefficient of V 3 is negative, the time rate 
of change of V 2 +fj 2 is negative and the conclusions of asymptotic stability of (1 1) follow. 

In order for the coefficient of V 3 to be negative, with the assumption that C D is positive, which is borne out in 
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This places limitations on the size of the elastic deformation and velocity. 

As has been observed earlier, the dynamics of q and £are decoupled from the dynamics of velocity (w, w,fj ) . 
Therefore, while studying the internal dynamics, it may be assumed that q and 6>are known functions and hence (11) 
becomes a three dimensional time-varying nonlinear system. 

Furthermore by design and physical limitation of the aircraft q ^ 0 and 6 — >0 cW as / — >qo , so the internal 
dynamics in (1 1) become 
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where F n F t e for i = x,z y rj have been substituted from (6). Mathematically, the residualized system (12) is called 
the limiting system of system (11) whose dynamics behavior determines that of (1 l) 5 and is discussed in detail in 
reference 6. 

Lemma 2: Define G as the right hand side of (12), F = (/T + F™ + f n } and M = diag(m'\m' ] ,m ~ l ) . Assume that 
the aerodynamic forces have a dissipative effect, that is 
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for all (u y w y rj) not equal to (0,0,0). Then the only possible orbits of the residualized system are equilibria. 


Proof: 

The proof is by contradiction and application of the Divergence Theorem. Assume that there exists a closed orbit C 
and let R define the interior region, then by the Divergence Theorem 

JjG ■ ndcr = JJJrf/VGrf Volume 

( R 

where tidcrxs, a vector element of surface area directed along the unit outer normal vector n to C. The inner product 
on the left side is zero by construction, while divG on the right is equal to - co ! + div(MF) and is negative by 
assumption. Hence, there exists a contradiction which proves the lemma. (This Lemma is an extension of the 
theorem in Hale (1980) chapter 2, Exercese 1.3.) 

As a corollary of the above dissipative condition lemma, the following stability result applies to the residualized 
system ( 12 ). 

Corollary: Assume that 

1. the residualized system (12) has a unique equilibrium for given T y 6 cmdy and GW*; 

2. the residualized aerodynamic force vector F satisfies the divergent condition (13). 
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Proof: By Lyapunov stability result of Lemma 1 any solution (z/(f), w(f), 77 (f)) is bounded. Thus, its limit set is 
either an equilibrium or a periodic orbit. But Lemma 2 excludes the periodic orbit case. Hence, the limit set of any 
solution consists of one point w 0 ,r] 0 ) which is the unique equilibrium of (12) by the assumption. This is 

equivalent of saying that ( u tj , w o , r } a ) is the global attractor of the system (12) 


Remark: For a given T, 0 cmd , and GW*; the residualized system (12) has a unique equilibrium if and only if the 
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have only one solution (a t) , VJ for given T y O cmdy and GW* results from residualizing (12) further to remove rj. 
Solution can be found by looking for an intersection of contour plots of each equation. 
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The culmination of lemmas and corollary combine into the main stability result. The stability result for the 
closed loop system by (10)-(1 1) is given in the following theorem: 
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Furthermore, any solution (u(t),w(t\q(t),0(t\ri(t),f}(t)) of the closed loop system (10)-(1 1) satisfies 
u(t)~*u (}9 w{t ) ^ w (? , ?(/) -> 0 , 0(t) -» 0 amh rj(t)^rj n , rj(t)-+ 0 

as t -» oo . In other words (m„, w„,O,0 c . mf/ ,7/„,O) is a global attractor of the closed loop system (10)-(1 1). 


Proof: A basic sketch of the proof is given here. A more rigorous and detailed version is presented elsewhere due to 
space limitation 6 . 

By the design of the feedback control the closed loop system is decoupled into internal dynamics (1 1) and controlled 
dynamics (10), which are stable. Furthermore, since q{t) 0 and 0{t) 0 cnui as t -* oo , internal dynamics system 
(1 1) is asymptotically autonomous and its limiting equation is given by the residualized system (12). We invoke the 
Markus theorem 7 to relate the stability of (1 1) to that of (12) which establishes that ( u () , w„ , rj ( ,) is the local attractor 
to the time varying system (11). Then we apply the Yoshizawa theorem 8 to show that it is a global attractor of (1 1). 

The first two assumptions are based on the physical characteristics of the vehicle under consideration. 
Furthermore, a variation of these assumptions that does not include the flexible vehicle dynamics, have been shown 
by Morton, et. al. to hold for a variety of fighter aircraft. Hence, the stability of the closed loop system depends on 
whether equations (14) have a unique solution (a o ,V 0 ). This in turn depends on the uniqueness of a for a given 
engine thrust T and the commanded attitude 0. 

In applying the stability results to the HSCT aircraft we take into consideration the structure of the aircraft 
model. The dynamic aeroelastic data is derived in a linear fashion and is then combined with nonlinear rigid body 
aerodynamic model 9 . The implications of this nonlinear/linear conglomeration on system shown in (1) is the 
removal of functional dependence of all q and rj terms on V. They all become linear functions of a. 
Define E {r) as a linear representation for the aerodynamic coefficients C (m) , then the following replacements are 
made in system (1): 

z n =^pSF 2 C, v (a,M,Gr*) =^pSV 2 C llj (a,M,GW*) 

M n =±pSV 2 cC u „(a,M,GW*) M„ = X -pSV 2 cC M 

E n =^ P SV 2 qC nn (M,GW*) E„ =^pSV 2 qC, in (M,GfV*) 

E u u + E w w+ E q q = ^pSV 2 (C v ,(M,GW*)u +C l?ll (M,G^*)w+C w (M,G^*>/) 

=± P SV 2 C^(M,GIV*) 

The aerodynamic functions for both nonlinear rigid body dynamics and dynamic aeroelastic contribution are shown 
in figures 2 and 3. Taking advantage of the model structure, assumption (1) of the stability theorem reduces to 

c„ (a) = C D (a) - C ns (a) (a) > 0 

and assumption (2) becomes 
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Figure 4 illustrates how the assumptions are satisfied for a sample flight condition chosen to produce the worse case 
dynamic aeroelastic interactions. In assumption 2, the coefficients of the all v\ and rj terms are several orders of 

magnitude smaller than Ef{a ) . By model construction m t? = 1 , thus making E™(a) a completely dominant term 
on the right hand side of the inequality. From a strictly mathematical perspective, one might argue that r| and 
especially 77 can become large enough to rival the dominance of E™ (a ) . If that were the case, then in the physical 

world the vehicle would have suffered a critical structural failure. For illustrative purposes, we assume 

77 = IOOO 77 =~ 120 ft! sec and show in figure 4 the difference between left and right hand sides of the inequality. 


Conclusions 


In this paper, we started to establish conditions under which the dynamic inversion control methodology can be 
guaranteed global stability when applied to a real world problem. The problem under consideration is an advanced 
high-speed, large flexible aircraft that requires the issue of flexibility to be to addressed as part of flight control 
design. The initial approach is to simplify the problem, address the pitch axis dynamics and control only one flight 
variable. The role played by flexible dynamics is immediately apparent from consideration of the internal dynamics 
of the system. Furthermore, these flexible dynamics play a role in establishing stability guarantees for the closed 
loop system. 

The results presented are the first to include flexible dynamics in stability analysis of dynamic inversion 
methodology. These form an initial basis to more complicated control problem formulation that includes an 
integrated structural/flight control that is essential for this class of vehicles. The results of this work will be 
presented in subsequent publication. 


Figures 





Figure 1. HSCT candidate vehicle configuration. 
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Figure 2. Nondimensional Aerodynamic Functions (a) Lift, Drag and Pitching Moment Coefficients and (b) Lift, 
Drag, and Pitching Moment Control Coefficients 
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Figure 3. Dimensional dynamic aeroelastic contributions to Lift, Drag, and Pitching Moment due to elastic 
deformation, 77 , and rate of change of the deformation, daJdt =tj * 
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Figure 4. Assumptions 1 and 2 
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